Magnetic field induced transitions in multiferroic TbMnOs probed by resonant and 

non-resonant X-ray diffraction 
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Multiferroic TbMnOs is investigated using x-ray diffraction in high magnetic fields. Measurements 
on first and second harmonic structural reflections due to modulations induced by the Mn and Tb 
magnetic order are presented as function of temperature and field oriented along the a and b- 
directions of the crystal. The relation to changes in ordering of the rare earth moments in applied 
field is discussed. Observations below T^ 6 without and with applied magnetic field point to a strong 
interaction of the rare earth order, the Mn moments and the lattice. Also, the incommensurate to 
commensurate transition of the wave vector at the critical fields is discussed with respect to the 
Tb and Mn magnetic order and a phase diagram on basis of these observations for magnetic fields 
H\\a and H\\b is presented. The observations point to a complicated and delicate magneto-elastic 
interaction as function of temperature and field. 
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PACS numbers: 75.47.Lx, 75.50.Ee, 77.80.Bh, 78.70.Ck 



I. INTRODUCTION 



Magneto-electric materials or multiferroics have stim- 
ulated much interest both scientifically and technologi- 
cally. Scientifically because the mechanism of coupling 
between ferroelectricity and magnetism is a fundamental 
part in understanding the properties of materials. Tech- 
nologically because of the flexibility of controlling states 
of a device with either electric or magnetic field or both. 
Although modern multiferroics operate at temperatures 
and fields prohibitive of direct applications they offer an 
exciting play ground in order to better understand such 
coupled behavior. 

Manganite perovskites i?Mn03 where R is a triva- 
lent rare earth ion have shown an extremely rich (H,T) 
phase diagram of ferroelectric and magnetic phases. 
Here the frustration of magnetism in i?Mn0 3 mangan- 
ite perovskites^ offers the means to break the chemical 
incompatibility between ferroelectricity and magnetism^ 2 - 
The tuning of the tolerance factor with appropriate size 
i?-ions suppresses Tjv for A-type ordering of ferromag- 
netic Mn layers that are stacked anti-parallel along the 
c— axis and leads to low temperature incommensurate 
(IC) spin ordering for R =Gd, Tb and Dyi For TbMn0 3 
in particular, Mn spins order below Tjv = 41 K to 
form a spin density wave (SDW) with propagation vec- 
tor T Mn = <5 M "b*, 5 Mn ~ 0.27, while below T s = 29 K 
the observation of a spontaneous polarization (P) along 
the c— axis coincides with the onset of a transverse spi- 
ral (cycloidal) ordering of Mn-spins£ At lower tempera- 
tures there is an additional transition below = 7 K 
in which Tb-spins order also incommcnsurately with 
r Tb = 5 Tb b*, 8 Tb = 0.42. Below T s , neutron diffrac- 
tion measurements show that the magnetic ordering is 
that of a transverse spiral with Mn spins rotating within 



the be— plane. Within the orthorhombic Pbnm crystal 
structure of TbMn0 3 (a = 5.3 A, b = 5.8 A, c = 7.4 A 
at 300 K), this ordering can be described as: 
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c x + m y ycos(T ■ r) + m z zsin(T ■ r) (1) 



where m s represent the magnitudes of the Mn moment 
along the principal crystallographic directions x, y, z 
(a, 6, c) and r is the position of the Mn ion. For zero field 
m x is very small or zero while below T s the orthogonal 
components along y and z result in a transverse spiral 
ordering of Mn-spins4 As discussed in Ref [H, from phe- 
nomenology it follows that the ferroelectric polarization 
is given by: 



JXem x m z [e x x r] 



(2) 



where %e is the dielectric susceptibility, 7 a coupling 
constant and e x is the spin rotation axis. Since the propa- 
gation vector is parallel to the 6-axis and Mn-spins rotate 
within the frc-plane around the a-axis (e K ||a), the above 
relation predicts the direction of the polarization to be 
along the c— axis as indeed is found expcrimentallyi^ Sim- 
ilarly the flop of the polarization from P\\c to P||c£ that 
occurs under magnetic field applied cither along the a— 
or b— axis, would be expected to result from a flop of the 
transverse spiral from the be— to the ab— planed 

The incommensurate magnetic ordering in multiferroic 
manganites is accompanied by lattice deformations that 
result in structural super-lattice reflections at 2tX The 
extinction condition in reciprocal space of these second 
harmonic reflections are the same as for their magnetic 
t = (5b* counterparts. Phenomcnologically the nature 
of these reflections is magnctostrictive and arises from 
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a quadratic magneto-elastic coupling! between an ampli- 
tude modulation of the magnetic moment and the lattice. 
However Jia et al. — have recently discussed the spin- 
lattice coupling in multiferroic manganites in terms of 
the magnetostrictive, orbital and Dzyaloshinskii-Moriya 
(DM) interactions (or Si x Sj spin-current terms^) on 
the basis of the electronic configurations of insulating 
manganites. From these results it follows that for a 
spin density wave (SDW), 2r lattice reflections are ex- 
pected to be observed, while for a circular transverse 
spiral [m y — m z ), 2r reflections arc suppressed.— If the 
transverse spiral becomes conical (m x ^ 0) it is possible 
to observe structural reflections also at t4 This means 
that first harmonic reflections can have a mixed mag- 
netic and lattice character. Indeed we have observed in 
magnetic fields the S Mn h* reflection in non-resonant x- 
ray diffraction experiments^ for {jLqH > 1 T. Therefore 
the investigation of incommensurate reflections in mul- 
tiferroic manganites using X-ray diffraction can provide 
valuable information on the type of magneto-elastic cou- 
pling that can be active in these multiferroics. 

More recently it was shown that in TbMnC>3 the mag- 
netic ordering of Tb- and Mn-spins are highly coupled 
below T s iii While for T > T s the magnetic wave vectors 
for Tb and Mn are locked so that r Tb = r Mn below 
it is found that r Tb and r lock-in to wave vectors 
whose magnitudes are rational fractions 3/7b* and 2/7b* 
respectively, while the wave vectors hold the relation- 
ship 3r Tb — T Mn = 1. This novel matching of wave vec- 
tors can be described within the frustrated Anisotropic- 
Next-Nearest-Neighbor-Ising (ANNNI) model coupled to 
a periodic external field produced by the Mn-spin or- 
der, as detailed in Ref. LLU. Within this model the 
r Tb = T Mn behavior is recovered while the r Tb = 3/7b* 
and r Mn = 2/7b* regime is stabilized by an optimal or- 
dering of 6 domain walls in the Tb spin-density wave, su- 
perimposed on the Mn-order. This model further shows 
that the ordering of r Tb = 3/7b* is energetically more 
favorable than the simple ft 1 1 with r = l/2b* found at 
low temperatures for DyMnC>3.— 

In this paper we report on x-ray diffraction measure- 
ments of the lattice deformation in TbMnC>3 as a func- 
tion of temperature and magnetic field applied parallel 
to the a— and b— axis. The paper is organized in the 
following way. In section II we present the experimental 
setups used for the present investigations. In section III 
we report on zero field measurements as a function of 
temperature for selected 2t reflections. In section IV, we 
present measurements with magnetic field applied along 
the a*-direction and in section V, with magnetic field ap- 
plied along the b*-direction. The phase diagram resulting 
from our data is discussed in section VI. 

In this paper we follow the convention of labeling in- 
commensurate magnetic reflections from Bertaut's repre- 
sentational theory.— Here the magnetic ordering of Mn- 
spins within space group Pbnm and t m "=0.27&* can be 
described by four irreducible representations (r) that 
consist of four modes (labeled as A, C, F and G)J^ 



The modes correspond to magnetic superlattice reflec- 
tions that occur in different Brillouin zones with extinc- 
tion conditions as follows; A for h + k=cvcn, Z=odd; G for 
h+k=odd with ^=odd; F for /i+fc=even and Z=even; and 
C for h + fc=odd and Z=even, where h, fc, I are Miller in- 
dices. Each mode here describes the Mn-spin polarization 
along one unique principle crystallographic axis.— It is 
shown that the Mn transverse spiral ordering in TbMnC>3 
is described by two irreducible representations T2 x T3 
where, T 2 = (C x ,F y ,A z ) and T 3 = (G x ,A y ,F z )£ The 
ordering of the Mn-spins is described by two A-modes, 
A y ,A z £ Using a combination of unpolarized and po- 
larized neutron diffraction and resonant magnetic X-ray 
scattering it has been shown that the remaining modes 
arise from the ordering of Tb-spins.— ^ In the following, 
the notation T Mn / Tb = (0, § Mn / Tb , 0) for first harmonic 
and 2r Mn ' Tb = (0, 2S Mn / Tb , 0) for second harmonic Mn 
and Tb incommensurate superlattice reflections is used, 
with wave vector r and wave number S. In the interme- 
diate and the commensurate phases, numerical values for 
wave numbers are used. 



II. EXPERIMENTAL 

The TbMn03 crystals used in our experiments were 
grown at the Hahn-Meitner-Institute in Berlin using the 
floating zone technique under Ar atmosphere. The crys- 
tals cut from the crystalline boule show an excellent crys- 
tal quality with a mosaic spread of 0.016° of the (0 2 0) re- 
flection (FWHM). In Fig. Q] we show measurements of the 
temperature dependence of heat capacity and magnetic 
susceptibility from a small single crystal cut from the 
same crystalline boule. The data shows three successive 
transitions with decreasing temperature at Tn = 41 K, 
T s = 28 K and = 7 K. The measurements and transi- 
tion temperatures are in good agreement with published 
measurements,— indicating a high quality crystal. 

The experiments were performed both at the beam- 
line BW5 at the Hamburger Synchrotronstrahlungslabor 
(HASYLAB) and at the beamlinc X21 at the National 
Synchrotron Light Source (NSLS) at Brookhavcn Na- 
tional Laboratory. 

At HASYLAB the experiment was conducted at a pho- 
ton energy of lOOkeV in horizontal scattering geometry. 
The sample was mounted in a Cryogenics superconduct- 
ing cryomagnet with horizontal field up to 10 T. The 
beam was monochromized by a (lll)-SiGe gradient crys- 
tal. A second SiGe gradient crystal was used as analyzer 
to suppress background. The single crystalline TbMn03 
sample had a size of 2x3x0.6 mm 3 and was measured 
in transmission geometry, which means the true bulk of 
the crystal is probed. The sample thickness in the di- 
rection of the beam was of the order of the absorption 
length at this high x-ray energy and was thus optimum 
for obtaining maximum scattering intensity. The crystal 
was aligned with the frc-plane in the horizontal diffraction 
plane in order to access (0, k, I) reflections. Using this 
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FIG. 1: (Color online) Temperature dependence of the spe- 
cific heat (C p ) and magnetic susceptibility \ of a small single 
crystal from the TbMn03 boule. For the magnetization mea- 
surements a field of 0.01 T was applied along the c— axis. The 
C P (T) and x(T) data indicate three successive phase transi- 
tions at 7 K, 27 K and 41 K corresponding to T™, T s , T N , 
respectively. 



setup, measurements were performed with field H\\b. 

At NSLS, the experiment was conducted in the hard 
x-ray regime at 9.5 keV as well as at the Tb absorp- 
tion edge with a photon energy of 8.252 keV. The sample 
was mounted in a 13 T Oxford cryomagnet with vertical 
magnetic field. With this setup, measurements with field 
H\\a were conducted, with the fee-plane oriented in the 
horizontal diffraction plane. 

A (0 2) graphite analyzer was used for background 
reduction for the non-resonant measurements. As detec- 
tor, a field insensitive Avalanche Photodiode (APD) was 
used. 

For measurements performed at the absorption edge 
polarization analysis was performed using the (0 6) 
reflection of the graphite analyzer. The scattering ge- 
ometry was 7r — a' with the analyzer at an angle of 
2#(oo6) — 84.6° implying a leakage of less then 1% from 
the 7T — 7r' channel. 

The sample used here was a crystal with polished c- 
surfacc and a size of about 3x3 mm 2 in order to obtain 
high diffracted intensities in Bragg geometry. 

All results presented in this study measured away from 
the resonances at the absorption edges are due to pure 
charge scattering. This can safely be assumed since the 
non-resonant magnetic scattering cross-section is about 
six orders of magnitude weaker than the charge scattering 
cross-section. Besides this, the magnetic signal of the 
observed reflections is further reduced due to the small 
magnetic form factor at the high Q- values we investigate 
here. 
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FIG. 2: (Color online) Temperature dependence at zero field 
and 9.5 keV photon energy of the (a) intensity and (b) k- 
component of the wave vector of the A- mode (0, 2+2S Mn , 
9) and F-mode (0, 2+25 Mn , 8) reflections. The respective 
temperatures Tjv, T s and are indicated by dashed lines. 
In addition data for the A-mode (0, 2S Mn , 3) reflection mea- 
sured at 100 keV are plotted. Due to a systematic offset the 
absolute wave numbers are not accurate. 



III. MEASUREMENTS AT ZERO FIELD 

In Fig. [2] we show the dependence of the intensity and 
wave vector with temperature of the A-mode (0, 2+25 Mn , 
9) and the F-mode (0, 2+2r5 M ", 8) reflections measured 
using a photon energy of 9.5 keV. From this data it is 
clear that the variation of intensity of the two reflec- 
tions with decreasing temperature is substantially dif- 
ferent and reflects the origin of these reflections. The 
A-mode reflects the lattice modulation that arises purely 
from the ordering of Mn-spins, while the F-mode reflects 
the lattice distortion that arises from the induced or- 
dering of Tb-spins with r Tb = r . For the A-mode 
reflection with decreasing temperature below we find 
a linear increase in its intensity up to T s while below 
T s its intensity decreases with further cooling down to 
T™. As Tb-spins order, the intensity of these reflections 
rapidly increases for T < T™. This is in sharp con- 
trast to the F-mode reflection where the intensity shows 
a smooth increase with decreasing temperature down to 
T^ b while a similar jump in intensity as for the A-mode 
reflections is observed below T™. In Fig. [2]we show also 
data measured from the A-mode (0, 2S Mn , 3) reflection 
using 100 keV x-rays. The behavior here is different to 
the 9.5 keV data. The intensity increases smoothly with 
temperature down to T^ h , below which it saturates. No 
jump in intensity nor a jump of the wave vector at 
is observed here. A possible scenario for the behavior of 
this reflection is presented at the end of this section. 
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FIG. 3: (Color online) Tb superlattice reflection as function 
of temperature, (a) shows the intensity of the (0, 2 + 2S Tb , 9) 
and the (0, S Tb , 3). (b) shows the variation of the propagation 
vector of the (0, 2 + 2S Tb , 9) superlattice reflection. In the 
inset, a scan over the first and second harmonic peak along 
the b*-direction is shown. 



In Fig. [3] we show detailed measurements of the (0, 
2 + 2S Tb , 9) and (0, 5 Tb , 3) reflections below = 
9 K that describe the lattice modulation associated with 



the change of the Tb-ordering from r 



Tb 



to 



37- 



Th 



Mr, 



=ljii The intensity of the 2r Tfc reflection 
shows a typical order parameters behavior with cooling 
below that follows the ordering of Tb-spins. Below 
6 K we can also measure the first harmonic (0, S Tb , 3) 
reflection (Fig. [3jb)) that is visible with non-resonant 
x-ray scattering showing that below there is a struc- 
tural component to the first harmonic magnetic reflec- 
tion. This is supported by the observation of a strong 
first harmonic Tb reflection shown in Fig. 2] 

The variation with temperature for the wave vectors 
2r Mn and 2r Tb is shown in Fig. [DJb) and [3£b) respec- 
tively. We find that the value of 28 Mn rapidly decreases 
on cooling through the SDW regime as noted earlier and 
shows a weaker temperature dependence below T s £■ How- 
ever at Tj^ b we note a substantial increase in the value 
of 2S Mn , that tracks closely the rapid changes in S Tb be- 
tween 9 and 7 K, while below 7 K the values of both in- 
commensurabilities remain relatively temperature invari- 
ant. Below T^ b we confirm the observation that the wave 
vectors for Mn- and Tb-spin order are coupled. Indeed 
the wave vectors of the (0, 2+2S Tb , 9) and (0, 2+25 Tb 7 
8) as well as the (0, 2+S Mn , 8) and (0, 2+S Mn , 9) in 
Fig. [4})-e are following the relation 3r Tb — r Mn = 1. Be- 
low 5 K the respective wave numbers approach the values 
of 5 Mn = 2/7 and S Tb = 3/7 within 0.002 accuracy, as 
was shown to be the case in Ref. UTI . 

We note here that we found surprising differences in 
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FIG. 4: (Color online) Temperature dependence of (a) inten- 
sities and (b-e) wave number of the first and second harmonic 
structural Mn- and Tb-superlattice reflections in the Tb or- 
dered phase T < T^ b at zero field. 



the change of the value of 5 Mn through the transition at 
T™. Through this transition Kenzelmann et al. (Ref. d) 
report a jump in S Mn of ~ 10~ 3 using single crystal neu- 
tron diffraction. This value is similar to the neutron mea- 
surements reported in Ref. 3 as well as X-ray diffraction 
measurements using 100 keV x-rays shown for the (0, 
2S Mn , 3) reflection in Fig. \Z[b). This is in contrast to 
the ~ 10 -2 change we find across using 9.5 keV x- 
rays in this work (Fig. [5Jb)) as well as in Ref. [ll|. In 
the case of the 9.5 keV experiments we probe few fim 
of the surface of the crystal while in the neutron and 
high energy x-ray experiments we probe the bulk of the 
sample. This one order of magnitude difference in the 
change of 5 Mn across suggests that strain effects at 
the crystal surface allows the lock-in of the wave vectors 
to values of rational fractions within AS ~ 10 -2 , while 
the unstrained bulk appears to modulate the values of 
the wave vectors below T™. 

We now turn our attention to the intensity variation 
with temperature of the A- and F-mode 2t m ™ reflec- 
tions shown in Fig. [2] The behavior of the A-mode (0, 
2+2S Mn , 9) reflection can be understood in terms of the 
changes in the magnetic ordering of Mn-spins. In the 
SDW regime between T s < T < Tjv the increase in in- 
tensity of the 2r M ™ A-mode arises from the quadratic 
magneto-elastic coupling in which the amplitude of the 
lattice displacements varies linearly with the increase in 
the size of ordered Mn moment The decrease in in- 
tensity below T s is ascribed to the development of a per- 
pendicular component that alters the collinear SDW to a 
transverse spiral^ For a perfect spiral with m y = m z , the 
2t m " reflection would be completely suppressed^ how- 
ever our data shows that at a temperature just above 
T^ b the 2t Mu reflection is still observed indicating that 
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the spiral remains elliptical (m y ^ m z ). Indeed this is 
confirmed by neutron diffraction by measuring directly 
the values of m y and m z to be 3.9 and 2.8 /is/A/n 
respectively^ At T^ b the intensity of the 2r Mn A-mode 
jumps dramatically. At first glance this rise in intensity 
may suggest changes in the magneto-elastic coupling of 
Mn- and Tb- sublattices. However, more likely the jump 
in intensity below arises from the matching of the 
Mn and Tb wave vectors as described above. Since we 
know that below we have significant structural con- 
tribution to the magnetic r Tb reflections it is most likely 
that here the r Tb intensity is superimposed on a 2r Mn re- 
flection below Tj^ b and thus providing dramatic changes 
in intensity. 

More puzzling is the difference in behavior of the A- 
mode reflections measured with 9.5 keV and 100 keV x- 
rays also shown in Fig. [2] where the 100 keV data show a 
linear variation in intensity with cooling from Tjv down to 
in contrast to the behavior observed for 9.5 keV pho- 
tons. Here again, we believe this has to do with surface 
versus bulk properties of the material, since also measure- 
ments performed with bulk sensitive neutron scattering 
show the behavior of the A-mode reflections we observe 
here with high-energy x-raysi&M 



IV. FIELD ORIENTATION H\\a 

The application of magnetic field along the a— axis re- 
sults in strong modulations of the magnetic order of Tb- 
and Mn-spins with the consequence of affecting signif- 
icantly the ferroelectric properties. The most notable 
change with field is the flop of the direction of the spon- 
taneous polarization from P\\c to P\\a. Phcnomcnologi- 
cally this is attributed to the flop of the Mn spin spiral 
from the be— plane to the ab— plane. This flop is accom- 
panied by a change of t Mu to a commensurate value 
of 1/4 at a critical field Hq of approximately 10 T at 
2 K In this section we focus predominantly on field in- 
duced magneto-elastic transitions that occur below Hq 
and give us an insight into both the magnetic ordering 
and magneto-elastic coupling in TbMnC>3. 



A. Phase region T% b < T < Tff n 

We measured the field dependence for H\\a for r and 2r 
reflections at T = 27, 15 and 9 K. In Fig.[5}r-b and Eb-d, 
we show the field dependence at T = 27 and 15 K, below 
T s , for the A-mode 2t m " reflection (0, 2 + 2S Mn , 9) and 
the F-mode (0, 2 + 2S Mn , 8) reflection for 15 K only. In 
this regime Tb- and Mn- spins are ordered with the same 
wave vector ( T Mn =T Tb ). The r M " and 2r M " reflections 
for the C- and G-modes are too weak to be observed at 
these two temperatures. 

The intensity of the 2r Mn A- and F-modes for T =27 
and 15 K is essentially constant with increasing field 
up to suggesting that field does not change cither 



(a) TTTTTTTT 



2.65 



2.45 



• •••••• 



• (0 2+28 Mn 9) H||a 
T(0 3-8 Mn 9) T=27K 



2.25 



<2 6 



H ^ 

M 




(C) T 



▼ ▼ ▼ 



• • • • • 



• (0 2+28 Mn 9) H||a 



A(0 2+28" 



T=15K 




02468 10 02468 10 
H H (Tesla) |i H (Tesla) 

FIG. 5: (Color online) Magnetic field dependence of (a) wave 
vector and (b) intensity of the second harmonic (0, 2 + 
2S Mn , 9) (circles) and first harmonic (0, 3 - S Mn , 9) (tri- 
angles) superlattice reflection at a sample temperature of 
T — 27 K for H\\a. Magnetic field dependence of (c) wave vec- 
tor and (d) intensities of the CM and IC superlattice reflection 
for magnetic field H\\a at a sample temperature T = 15 K. 
At (j,qHq = 9 T the reflections lock in at CM positions. The 
horizontal dashed lines show the half and quarter integer po- 
sitions of the CM first and second harmonic reflections above 
H a c . 



the ellipticity of the Mn spin-spiral or the magneto- 
elastic coupling associated with it. Surprisingly above 
2 T for both 27 and 15 K we observe the T Mn G-mode 
(0, 3— 5 Mn , 9) reflection. The wave number of this reflec- 
tion is <5 M "=0.275 at /iqH = 2 T, in agreement with the 
incommensurability 2<5 M "=0.55 of the 2r M " reflection. 
The intensity of this reflection shows a linear increase 
with field above 2 T and saturates at approximately 6 T 
(Fig. [5jj and d). Similar measurements as a function of 
field at 9 K just above show the very same behavior 
as described above for the same reflections. Additional 
scans (not shown) demonstrate that only r Mn G- and 
C-modes are observed above 2 T, while the 2r Mn A- and 
F-modes are always measurable irrespective of field. The 
observation of t m ™ reflections in fields larger than 2 T 
may have a number of possible reasons. The emergence of 
a first harmonic reflection may be thought at first glance 
as trivial as the coupling of the lattice to H becomes lin- 
ear. However Jia et al. suggest that a linear behavior to 
the spin-lattice coupling may occur in the case when the 
transverse spiral becomes conical i.e. m x ^ 0£ If this 
was indeed true then we would expect to observe the 
r Mn reflection of an A-mode as such reflection describe 
the Mn spiral ordering. However according to our obser- 
vations the first harmonic reflections in field arc found for 
G- and C-modes suggesting possible changes to the Tb- 
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spin ordering with magnetic field. Indeed magnetization 
measurements in Ref. [15J suggest the ferromagnetic align- 
ment of Tb-spins with H\\a for T=15 and 9 K. Therefore 
the emergence of r Mn reflections with field most likely 
is associated with changes in the Tb-ordering as opposed 
to changes in the magneto-elastic coupling. 



B. Phase region T < 

We now turn our attention to the behavior of the in- 
commensurate reflections as a function of field below 
and remind the reader that in this regime at zero field 
the Tb and Mn magnetic ordering is coupled so that 
3r Th _ T Mn = i. As shown Fig. [6] we find r Tb and the 
A-modc r Mn for n H = 1 T and T < T™. Increasing 
magnetic field results in the strong attenuation of the 
2r Tb (0, 2+2S Tb , 9) reflection and its complete suppres- 
sion above noH > 2 T, consistent with the FM ordering 
of Tb-spins as indicated by neutron diffraction*^ 

For intermediate field values (0 < H < 2 T) the 2r Tb 
reflection vanishes from the 25 Tb = 0.86 position and 
seems to shift to 0.90 (Fig. [5]). On the other hand reso- 
nant scattering from the Tb L 2 absorption edge indicates 
a shift of the first harmonic r Tb reflection to 0.364 indi- 
cating a shift of the 2r Tb reflection to the 0.727 position, 
where also a peak appears at wl T as will be shown later. 
This reflection nevertheless stays also above a field of 2 T, 
which indicates that we deal here with the r Mn reflection 
observed also at temperatures T > T™. The shift ob- 
served at 1 T is coupled to a discontinuous change in 
the 25 Mn reflection from 0.572 to 0.56 also for the same 
field values as shown in the inset of Fig. E^b). This again 
is due to the superposition of the r Tb intensity on the 
2r Mn reflection at low fields and a pure 

2 T Mn re fi ec tion 

at higher fields. 

To investigate the transition to this intermediate phase 
we performed resonant X-ray diffraction measurements 
with and without polarization analysis in the n — a' 
polarization channel as a function of H\\a, by tuning 
the photon energy to the Tb L2-edge at 8.252 keV. 
At this photon energy, the background was significantly 
higher than for the experiments at 9.5 keV and prohib- 
ited the observation of any non-resonant scattering sig- 
nal. However, a strong resonant magnetic signal was 
observed at T = 4 K for the first harmonic Tb reflec- 
tions (0, k±d Tb , I) with I = odd ( Tb magnetic A- and 
C- modes). The measurements were conducted with an- 
alyzer used to reduce background (PG(002)) and as po- 
larization analyzer (PG(006)) for the tt — a 1 channel. In 
the inset of Fig. [7l energy scans with and without po- 
larization analysis are shown. The small intensity in the 
7r — er'-channel is due to leakage since the analyzer angle 
is not exactly at 29 = 90° but at 26» (006) = 84.6°. All in- 
tensity is thus scattered in the tt — tt' channel. According 
to the resonant magnetic scattering cross section, this is 
the case if the magnetic Tb moment is aligned along the 
a*-direction perpendicular to the scattering plane, con- 
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FIG. 6: (Color online) Field dependence at T = 4 K of the 
terbium and manganese superlattice reflections with magnetic 
field H\\a. (a) peak positions (closed symbols) and (b) the 
corresponding intensities (open symbols) are shown. Open 
rectangles show the (0, 0, 4) main Bragg intensity divided by 
a factor 100. The inset in (b) shows the peak profiles of the 
(0, 3-8 Tb , 9) and (0, 2+25 Tb , 9) reflections for different fields. 

sistent with results in Ref. [l^. 

In Fig. scans over the (0, 5 Tb ', 5) and (0, l-5 Tb , 
5) positions with 8 Tb = 0.43 are shown as function of 
field. As the field is increased, the intensity of these 
reflections decreases and the wave number S Tb slightly 
shifts as observed in the non-resonant experiment. At 
HqH = 0.9 T, this resonant reflection shifts discontinu- 
ously from S Tb = 0.428(2) (fi H = T) to S Tb = 0.363(1) 
(HqH = 1 T). While the wave number of this reflection 
is invariant in field, its intensity is rapidly reduced with 
increasing field and vanished above H\\a = 2 T, when 
Tb-spins show a FM ordering ordering^ Our resonant 
and non-resonant measurements demonstrate that this 
intermediate phase is stable between H\\a = 0.9 and 2 T. 

The discontinuous transition we find here at 1 T sug- 
gests a change in the coupling of the Mn and Tb ordering. 
In this new regime the values of 8 Tb and S Mn are close to 
the rational fractions 4/11 and 3/11, respectively, while 
the wave vectors hold the relationship of 2r Tb + r = 1 
with 0.001 accuracy. Modeling of the Mn- and Tb-spins 
using an ANNNI model it was shown that such a state is 
stabilized as the Tb-SDW acquires a homogeneous com- 
ponent via \H an( i couples to the quartic term of the 
Landau expansion.— Interestingly the change in the cou- 
pling between Mn and Tb magnetic ordering character- 
izes a region in which the value of P\\c increases by ~30% 
up to H\\a = 2 T. 

For H\\a > 2 T the FM alignment of Tb-spins results 
in the melting of the incommensurate Tb-orderingi 4 i 10 i 15 
This melting is associated with a decrease of P\\c to val- 
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FIG. 7: (Color online) K-scans of the (0, S Tb , 5) and (0, 1-<5 T \ 
5) reflection as function of field H\\a at a sample temperature 
T = 4 K without polarization analysis. The inset shows an 
energy scan over the Tb L2 absorption edge with and with- 
out polarization analysis in the it — a' channel. The total 
signal only shows the resonant intensity. The fluorescence 
background is subtracted using an energy scan performed off 
the reflection position. 
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FIG. 8: (Color online) Incommensurability S of the super- 
lattice reflections (0, 2S Mn , 4) and (0, l-2S Mn , 4) below H h c 
(closed circles) and (0, 0.25, 4),(0, 0.5, 4) and (0, 0.75, 4) 
above Hq (open circles), measured at (a) T = 2 K, (b) 
T — 10 K and (c) T = 20 K as function of magnetic field 
H\\b. At T = 2 K in addition (0, S Tb , 5) and (0, 2S Tb , 5) 
reflections are shown below and above (open triangles). 



ues 30% smaller than the maximum polarization at zero 
field. However r Mn and 2t " reflections are observed 
above 2 T and discontinuously shift to values of 1/4 and 
1/2 respectively at H c (not shown). 



V. FIELD ORIENTATION H\\b 

Whereas the H\\a measurements were performed with 
photon energies of 9.5 keV that allows to probe only the 
first few fim of the crystal, the measurements with H\\b 
are performed with high photon energies of 100 keV with 
which the true bulk of the crystal is probed. For the H\\b 
geometry the sequence of magneto-structural transitions 
towards the polarization flop transition is significantly 
different than what we have found for the H\\a configu- 
ration. Although we find that the polarization flop again 
corresponds to a first order transition to a commensu- 
rate phase with <5 Mn =l/4, the sequence of transitions 
associated with Tb-ordering differs from the H\\a case. 
A summary of our results with respect to the wave num- 
bers S Mn and S Tb is shown in Fig. [HJi-c for T = 2, 10 
and 20 K respectively. The transition to 5 Mn = 1/4 is 
clearly evident from the data while Hq decreases with 
decreasing temperature to values that agree well with 
the observation of the flop in polarization from P\\c to 
P\\a&&21 
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FIG. 9: (Color online) (a) Intensities as function of field of 
the IC (0, 5 Mn , 3) and (0, 2S Mn , 4) reflections below H b c 
and the CM (0, 0.25, 4), (0, 0.5, 4) and (0, 0.75, 4) above 
H b c . (b) The respective wave vectors (0 k 4) for the intensities 
in (a) are shown. 



A. Phase region Tjj b <T < T^f n 

In the region ofT^ b <T<T N , we find only reflections 
with 2t m " as expected from previous measurements and 
indicative of quadratic magneto-elastic coupling. How- 
ever above H% > 1 T we find also first harmonic reflec- 
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tion such as the F-modc (0, S Mn , 4). This reflection that 
arises from the ordering of Tb-spins with r Tb = r Mn 
shows a linear increase in intensity with fields up to Hq 
indicating changes in the magnetic ordering of Tb-spins. 
In contrast its second harmonic counterpart (0, 2<5 M ™, 4) 
remains constant in intensity up to H c (Fig. [9^) • This 
linear increase over the whole field range is different from 
the behavior found for H \\ a (section IIV|) , where the in- 
tensity of these reflections increases initially but then de- 
creases again towards Hq. 

As we approach the critical field H h c the value of 2r Mn 
begins to move towards 1/2, while at H h c we find a first 
order phase transition where the intensity of the incom- 
mensurate reflection decreases with field and vanishes as 
the commensurate reflection appears (Fig. [9)3). With de- 
creasing temperature the value of H c also decreases as 
shown in Fig. [5h.-c. These changes to r Mn and 2r Mn F- 
mode reflections have to be attributed to changes in the 
magnetic structure associated with Tb-spins. We spec- 
ulate that at this high temperature the application of 
magnetic field may result in the polarization of the Tb- 
spins along the 6— axis. For this field direction there is no 
evidence of ferromagnetic ordering of Tb-spins as found 
for H\\a£&±2. 

Detailed measurements of the transition to the com- 
mensurate phase were made at 10 K by tracking the lat- 
tice modulation associated with the F-modc structural 
peak (0, 4 — 25 Mn , 0). These measurements, shown in 
Fig. QJJ were taken by cooling the sample to 10 K in zero 
field and applying field up to H \\a = 6.5 T and then down 
to T while the data were measured. These data show a 
hysteresis in terms of field and wave number 25 Mn that 
tracks the hysteresis in the polarization flop from P\\c to 
P\\a&m. In terms of magnetic field the hysteresis is as 
small as ~0.2 T. However in terms of wave number the 
effect is more pronounced. For increasing field the wave 
number remains at a value of 2S Mn = 0.557 almost up 
to H c , before it locks to the commensurate value of 1/2. 
For decreasing field, a rapid increase in 2S Mn just below 
Hq and a significantly lower value of 2S Mn = 0.543 at 
HoH = T is observed. Interestingly, these lower val- 
ues suggest S Mn = 0.272 which is close to the rational 
fraction 3/11 as found for the case of H\\a and 9.5 keV 
photons below T™. In the inset of Fig. [TD1 scans over the 
superlattice reflections for increasing field around H b c are 
displayed. It shows the coexistence of the incommensu- 
rate and commensurate phase right at the transition at 
HoH c = 5.5 T. The intensities of these reflection in the 
incommensurate phase for increasing and decreasing field 
show exactly the same values while only their positions 
show a hysteresis. 



B. Tb-ordered phase region T < T" 

The application of magnetic field at low temperature 
for the H\\b configuration results in a significantly dif- 
ferent behavior of r Tb and r Mn compared to the H\\a 
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FIG. 10: (Color online) Hysteresis of the superlattice reflec- 
tions measured at T = 10 K. Solid circles show the data 
measured with increasing field, open circles the data mea- 
sured with decreasing field. The inset shows the superlattice 
reflections at the phase transition for increasing field. 



configuration. Here the sample was cooled to 2 K in zero 
field and then data were measured with increasing field 
up to fj,oH\\a = 10 T, (Fig. (5£i). Between an applied field 
of 1 and 2 T, a phase transition is observed in the Tb 
sublattice with a shift of the peak positions from the in- 
commensurate value T Tb = 0.436* to the commensurate 
value T Tb = 1/36* and accordingly from 2r Tb = 0.866* 
to 2r Tb = 2/36*. Although this transition is only related 
to the Tb order, the intensities of the Mn reflections are 
influenced by this transition, whereas the wave vectors 
are not affected. As shown in Fig. [Til the intensities 
of the (0, 0.55, 3) (A-mode) and the (0, 0.55, 4) (F- 
mode) reflections show a distinct behavior. While the 
(0, 0.55, 3) intensity decreases at H™, the (0, 0.55, 4) 
intensity increases. 

Cooling the sample in an applied field [j,qH = 3 T 
we find in the temperature dependence that the transi- 
tion to the incommensurate phase with r Th = 0.436* is 
suppressed at = 4.5 K with decreasing temperature 
in preference for the commensurate r Tb = 1/36* phase 
(Fig. HJ). For the F-mode reflections (0, S Mn , 4) and 
(0, 2r5 Mn , 4) we find no changes in the value of S Mn but 
again a clear change in the intensity is observed at T™. 

In this field-cooled measurement below T = 6 K we 
find two previously unreported modulations at (0, 0.298, 
5), (0, 0.6, 5). The nature of these reflection is not clear 
at this point, however their field dependence indicates 
that they are associated with Tb spin ordering. We find 
that increasing field and increasing temperature result 
in a considerable decrease in their intensity (Fig. 113b - 
b). The reflections disappear at the transition to the 
<5 T6 =l/3 phase. Since these reflections only occur at the 
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FIG. 11: (Color online) Intensities of the (0 0.55 3), (0 0.55 4), 
and (0 0.66 5) reflections below H b c and (0 0.25 4) and (0 0.5 4) 
above Hq as function of magnetic field H\\b measured at a 
sample temperature of T — 4 K. 
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FIG. 13: (Color online) k-scans over (a) (0, 5, 5) and (b) 
(0, 26, 5) positions as function of magnetic field H\\b and 
temperature below ■ 



C. High-Field Commensurate Phase 




FIG. 12: (Color online) (a) Temperature dependence of the 
intensities of the Mn and Tb first and second harmonic su- 
perlattice reflections measured at magnetic field fioH = 3 T. 
(b) shows the corresponding k-components of the wave vector. 
The measurement is performed by going from the intermedi- 
ate Tb ordered phase through into the purely Mn ordered 
phase. 



high-energy x-ray setup they can be regarded as a bulk 
property related to Tb magnetic order below T 1 



Tb 
N ■ 



We now turn our attention to the commensurate phase 
above H h c . The scattering geometry that allows us to 
measure the commensurate phase makes only F- and C- 
modes accessible (Fig. [9]). At T = 10 K the transition 
to a commensurate phase occurs at ^.oH^ = 9 T. Above 
the F- and C- modes exhibit a linear increase of 
their intensity with field for the first harmonic reflec- 
tions, whereas the intensity of the second harmonic peaks 
remains constant with increasing field. This behavior 
follows on from the lower field behavior in the incom- 
mensurate phase. As the F- and C- modes arise from 
Tb-ordcring it only allows us to comment on the value of 
r Mn and not to changes associated with the flop of the 
Mn spin-spiral. Cooling the sample to T = 2 K and ap- 
plying field we find the transition to the commensurate 
phase at HoH^ = 5 T. Here reflections associated with 
the T Tb — 1/36* phase vanish at H h c and only commen- 
surate reflections with S Mn = 1/4 and 2r Mn = I /2b* are 
observed (Fig. [8]). The integrated intensities of (0, 0.25, 
4), (0, 0.5, 4) and (0, 0.75, 4) are plotted as function of 
H\\b in Fig. [TTJ Intensities of the 2t m ™ reflection are a 
factor of 30 weaker than the first harmonics. While the 
intensities of the (0, 0.25, 4) and (0, 0.75, 4) reflections 
increase linearly with field, the intensity of the (0, 0.5, 
4) remains constant up to maximum field. The Tb (0, 
0.66, 5) reflection also shows constant intensity from its 
onset at ^lqH = 2 T up to H h c . At H h c a little intensity 
is still present and coexists with the reflections of the 
commensurate phase. 

Following the F- and C-mode reflections while cooling 
the sample in a field /lo-ff = 6 T we find that r Mra first 
takes an incommensurate value which decreases with low- 
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FIG. 14: (Color online) (a) Intensity behavior of CM and IC 
superlattice reflections measured at H = 6 T below (CM) and 
above (IC) the critical temperature, (b) Corresponding wave 
number 5 as function of temperature for the first harmonic 
F-mode reflection, (c) k-scans in the region of coexistence of 
the CM and IC phase at the critical temperature. 



ering the temperature until T =T S (see Fig. ITIb.b). Be- 
low T s the value of T Mn remains relatively constant until 
approximately 18 K where it rapidly decreases and locks 
into a value of T Mn = 1 /4b* below 12 K. In Fig. [Ht, 
k-scans in the region of coexistence of the commensu- 
rate and incommensurate phase at the critical temper- 
ature are shown. The variation in the intensity of the 
commensurate reflections is quite different between the 
different reflections we examined. Whereas the weak sec- 
ond harmonic reflection is constant in intensity there is 
a strong variation for the first harmonic reflections, indi- 
cating changes in the Tb order are obviously still present 
even in the high field commensurate phase. 



VI. PHASE DIAGRAM AND DISCUSSION 

From the measurements we have described in this pa- 
per we can construct the H — T phase diagrams for H\\a 
and if || 6 (Fig- 115b and b) which summarize the various 
phases we have found and their correlation to the ferro- 
electric polarization of TbMnC>3 reported in Ref. [H. 

In zero field we find purely structural 2t m " reflections 
which arc interpreted to arise from quadratic magneto- 
elastic coupling. In the SDW regime T N > T >T S the 
intensity of the A-mode 2r Mn reflections rises with cool- 
ing as expected for quadratic magneto-elastic coupling 
and below T s it decreases as the Mn magnetic order be- 
comes spiral. While ideally 2r Mn reflections should dis- 
appear for a circular spiral^ their continual observation 
for T > T s reflects the ellipticity of the Mn spiral order- 
ing. 
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FIG. 15: (Color online) Phase diagram of TbMnOa for fields 
(a) H\\a and (b) H\\b. In (a) the light shaded area represents 
the intermediate phase region and the dark shaded area the 
Tb AF ordered state found for field H\\a. In (b) the green line 
represents as found for magnetic field H\\b. The green 
shaded area shows the intermediate CM Tb ordered phase for 
H\\b. The IC to CM transition for the manganese order is 
shown for H\\a and H\\b. The vertical dashed line defines the 
onset of the first harmonic reflections. 



In the regime T s > T > T Tb Mn- and and Tb-spin or- 
dering is coupled and T Mn =T Tb ±L Below T™, Tb-spins 
order with a different wave vector but remain harmoni- 
cally coupled to the Mn spiral ordering as their wave vec- 
tors assume values that obey the relation 3r Tb — r Mn = 1. 
This is confirmed in this work to an accuracy of 0.003 and 
we depict the phase regions of this coupled ordering in 
Fig. [T3i and b. 

Application of magnetic field in this harmonically cou- 
pled state results in different behavior depending on the 
direction of the field. When a relatively small magnetic 
field is applied along the a— axis below we find an 
intermediate state that is characterized by the appear- 
ance of a Tb reflection at 0.9 (Figs. SJd and [6^,), a jump 
of the wave vector r Mn to smaller values and a shift of 
T Tb to 0.366* observed by resonant magnetic scattering 



(Fig. [7]). The observed shifts in wave vector are com- 
pletely in agreement with the relationship 2r Th +T A/ " = 1 
(Sections HI] and llVB|) found in Ref. QJ. This region is 
marked as a light shaded area in Fig. 115a and appears 
both in the temperature dependence at zero field (Fig. 2]) 
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between Tjf = 7 K and ~ 9 K as well as in the field de- 
pendence at T = 4 K between 1 and 2 T (Fig. [BJ. For 
[XoH\\a > 2 T Tb spins become ferromagnetically polar- 
ized. This intermediate phase may be visible only in the 
surface near region probed in this field configuration and 
is likely to be absent if the true bulk of the crystal is 
probed. 

It is worth noting at this point that the intermediate 
value of T Tb below for both H \\a and H\\b configu- 
rations is in fact quite similar; r Tb =0.366* in the former 
case and T Tfc =0.336* for the latter. The difference in 
these two ordcrings of Tb-spins is that for H \\a r Mn is 
harmonically coupled to the Tb-ordering while for H\\b 
it is not. This subtle but significant difference we believe 
reflects the anisotropy of the Tb-spin ordering as well as 
the delicate balance of the coupling between Mn and Tb 
spins. 

When field is applied along the 6-axis below T^ b the 
harmonically coupled regime disappears and instead we 
find a phase transition of the Tb moments at about 
HoH = 1.25 T from the incommensurate to a commen- 
surate order with r Tb = 1/36* and the 2nd harmonic 
2r Th = 2/36* (sec Fig.HSb)- This intermediate Tb phase 
does not affect the wave vector r Mn , which only becomes 
commensurate at noH^, ~ 5 T. Nevertheless, intensi- 
ties are affected which emphasizes again the effect of the 
Tb magnetic order on the structure. Similar to H\\a, a 
dashed vertical line in Fig. [T5b shows the field of onset for 
the first order satellites which arc induced by the polar- 
ization of the Tb magnetic moment by the magnetic field. 
As discussed before (Section llVBj) in the Tb-ordered re- 
gion T < and the intermediate region (^H > 0), 
a clear assignment of Mn or Tb reflections is no longer 
possible. 

At higher field values for fj,oH\\a > 2 T and below T s 
only incommensurate Mn reflections with wave vectors 
r Mn and 2t m ™ are observed, whereas the Tb reflections 
disappear. Since first harmonic reflections in field are 
found for G- and C-modes only and not for the Mn-order 
related A-type reflections, this suggests that the vertical 
dashed line in Fig. [T5r marks the onset of the FM order of 
the Tb moments for applied fields H\\a. This also corrob- 
orates neutron diffraction measurements which find the 
sharp decrease of F- and C-modes in this field range and 
the sharp increase of ferromagnetic Bragg reflections. 10 
This FM ordering of Tb-spins is relatively stable up to 
fields H\\a > 8 T. 

For both field configurations H\\a and H \\b we find that 
the flop of the polarization is associated with a first or- 
der transition to an commensur ate T Mn = 1/46* phase. 
As the direction of the polarization changes, this flop is 
taken to reflect the flop of the spiral plane from the 6c— 
to the ab— plane. The jump to a commensurate wave 
vector at the spiral flop transition is not surprising given 
that 5 Mn is in the proximity of a commensurate value 
in zero field and therefore it is likely an energetically fa- 
vorable state can be reached by a relative small shift to 
S Mn = 1/4. This is in contrast to DyMn0 3 were the S Mn 



is not close to a commensurate value (r Mn — 0.386*) 
and at the polarization flop transition the wave vector 
does not shift at all^ In TbMnOs above H^, variations 
in Tb order seem still to be present as can be deduced 
from the strong variation of the first harmonic Mn reflec- 
tions as function of temperature at a field of 6 T shown 
in Fig. [T3J Here the transition of Tb ordering into the 
r Tb = 1/36* phase seems to extend into the commen- 
surate phase above Hq, as is shown by the horizontal 
point-dashed line in Fig. [15b . 

In more general terms the flop of the polarization from 
P\\c to P\\a can be understood by the magnetic degrees 
of freedom of the Mn transverse spiral. There are three 
excitations of the spiral. One is the phase of the spiral 
polarized within the 6c-plane. Two other components are 
polarized along the a— axis in zero fields One of these 
components results in a flop of the 6c-spiral by a rotation 
around the 6— axis while the other results in a twist of the 
spiral by a rotation around the c— axis. All measurements 
thus far are consistent with the assumption that field 
applied either along the a— or the 6— axis couples to the 
former mode and produces the flop of the spiral from the 
6c— to the a6— plane, while field applied along the c— axis 
melts the spiral ordering^ - 

As this picture can describe the main flop of the polar- 
ization as well as the magnetic excitations, this work has 
shown that the magnetic interactions between Mn and 
Tb spins remain significant and may be considered as a 
perturbation upon this model of the flop of the polariza- 
tion and the spiral. The strength of he Jj\/ n -i?, interaction 
can be gauged by noting that the flopping fields for the 
i?=Dy are lower than in the present case for TbMnO,-}. 
This can be understood as a weaker Jmu-r for the case 
of R=Dy. Indeed in DyMnOs at low temperature Dy- 
spins become completely uncoupled from Mn-spins and 
Dy orders with wave vector T Dy = l/2b*, a behavior sig- 
nificantly different from the harmonic coupling we find in 
TbMn0 3 below T^ b ±L 



VII. CONCLUSIONS 

In this paper we present an extensive study of the be- 
haviour of structural and magnetic superlattice reflec- 
tions of TbMnOs as function of temperature and applied 
magnetic field. Details in the H — T phase diagram are 
revealed by small changes in the magnitude of wave vec- 
tor and intensities. These findings help us explain the 
changes in spontaneous polarization as function of tem- 
perature and field. The subtle difference in the Tb order 
wave vector observed at low temperatures for non-zero 
fields H || a and H\\b reflects the anisotropy of the Tb- 
spin order and the delicate balance between Mn and Tb 
spins. This again highlights the significance of interac- 
tions between Mn and Tb spins for the description of the 
polarization flop at the critical fields. Differences in the 
magnetic wave vector in the crystal bulk and the surface 
near region, as is observed by high-energy (100 keV) and 
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hard x-ray investigations, show the influence of surface 
effects on magnetic polarization behaviour in multiferroic 
compounds. 
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